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Abstract 

We have determined the physical ground state properties of the compounds CeRuPO and 
CeOsPO by means of magnetic susceptibility x(T), specific heat C(T), electrical resistivity p(T), 
and thermopower S(T) measurements, x(T) reveals a trivalent 4/ 1 cerium state in both com- 
pounds. For CeRuPO a pronounced decrease of p(T) below 50 K indicates the onset of coherent 
Kondo scattering, which is confirmed by enhanced S(T). The temperature and magnetic field 
dependence of x(T) and C(T) evidence ferromagnetic (FM) order at Tq = 15 K. Thus, CeRuPO 
seems to be one of the rare example of a FM Kondo lattice. In contrast, CeOsPO shows antifer- 
romagnetic order at Tn = 4.4 K despite only minor changes in lattice parameters and electronic 
configuration. Additional 31 P NMR results support these scenarios. LSDA+C/ calculations evi- 
dence a quasi two dimensional electronic band structure, reflecting a strong covalent bonding within 
the CeO and RuP layers and a weak ionic like bonding between the layers. 
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I. INTRODUCTION 



Intermetallic Kondo lattice systems have attracted considerable attention in the last 
decades. In these systems, two different interactions determine the ground state properties, 
which can be understood in a simple way in terms of the Doniach picture.- The formation 
of a nonmagnetic Kondo singlet competes with the formation of a magnetic ordered ground 
state;- in the case of dominating Kondo interactions, the magnetic moments of the 4/ elec- 
trons are screened by the conduction electrons, and a nonmagnetic ground state develops. 
If the Ruderman-Kittel-Kasuya-Yosida (RKKY) exchange interaction exceeds the Kondo 
screening, long range magnetic order occurs. While many Ce-based Kondo lattices show 
antiferromagnetic (AFM) ground states, only very few systems are known with ferromag- 
netic (FM) order and pronounced Kondo effects. Therefore, the behavior at a FM critical 
point in a Kondo lattice is presently not settled. 3 From the few FM Kondo systems known, 
the best studied examples are YbNiSn, 4 CeSi2-a;,- and the compound-series CePdi_ x Rhj;.- 
However, in the former compound pressure stabilizes the magnetic state, preventing a study 
of the critical region. In the latter, disorder is not negligible and may hide the intrinsic 
physical properties. On the way to look for new Ce-based Kondo lattice systems close to a 
FM quantum phase transition, the CeTPO (T = transition metal) compound series attracted 
our interest, because of the rather unusual crystal structure with alternating layers of TP4 
and OCe4 tetrahedra, causing Ce-Ce inter-atomic distances to be in the range where FM 
order can be observed/ This compound series, crystallizing with the tetragonal ZrCuSiAs 
type structure, was discovered by the Jeitschko group;-^ but no physical properties have 
been reported for CeRuPO. The compound CeOsPO has not been synthesized before. Re- 
cently, LaFePO was found to be superconducting, supporting the interest in this class of 
material.— *ii 

In this contribution, we will present the low temperature physical properties of CeRuPO 
and CeOsPO. The article is organized as follows: In section Hi Al we will discuss the details 
of the sample preparation and the crystal structure. The physical properties of CeRuPO 
are presented in detail in section III Bl and compared with the results of CeOsPO in section 
III CI We present and discuss 31 P NMR results for both compounds in section III Dl Finally, 
in section III II we show band structure calculations of both compounds using an LSDA+[7 
approach and analyze the bonding properties and the magnetic exchange. 
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II. EXPERIMENTAL 



A. Preparation and Crystal structure 

The polycrystalline samples were prepared using a Sri-flux method in evacuated quartz 
tubes.— 2 For CeRuPO the preparation was optimized with Ru02 as the source for the oxygen. 
We found the best growing conditions for an initial molar ratio of 8:6:4:82 (Ce:P:Ru02:Sn). 
This preparation route was also used to obtain the LaRuPO compound, serving as the 
nonmagnetic reference. CeOsPO was prepared with SnC>2 instead of the hazardous osmi- 
umtetraoxide. For this compound we found the best results for the initial molar ratio of 
4:3:2:2:89 (Ce:Os:P:Sn02:Sn). The heating schedule was similar for all preparation cycles. 
The starting mixture was heated with a rate of 100°C per hour to 1000°C and kept at this 
temperature for 5 hours; subsequently, the melt was cooled down to 500°C by moving the 
quartz tube out of the furnace with an average cooling rate of 5°C per hour. After the 
reaction the excess Sn was dissolved in dilute HC1 which attacks the CeTPO at a lower rate. 
Small single crystals (< 50 jum) of CeTPO are attached to each other forming flake-like 
polycrystals which were ground and pelletized for further physical characterization. Several 
samples were investigated with energy dispersive X-ray analysis. These measurements reveal 
a stoichiometric Ce:Ru:P content and confirmed the presence of oxygen in the compounds; 
furthermore, no foreign phases could be detected. In addition, carrier gas-hot extraction 
(LECO, TCH 600) was used to determine the oxygen content xo of CeRuPO more pre- 
cisely. The result of iq = (25.0 ± 0.7) at% indicates a stoichiometric oxygen occupancy. 

Several powder X-ray diffraction patterns recorded on a Stoe diflractometer in transmis- 
sion mode using a monochromated Cu-K Q radiation (A = 1.5406 A) confirmed the formation 
of single phase CeTPO. The lattice parameters refined by simple least square fitting for 
CeRuPO (a = 4.028(1) A, c = 8.256(2) A) and LaRuPO (a = 4.048(2) A, c = 8.410(4) A) 
were found to be in good agreement with the reported single crystal data.- For CeOsPO, the 
first refinement was performed using an internal silicon standard yielding a = 4.031(1) A, 
c = 8.286(3) A, which was later confirmed in various measurements from different batches. 
The tetragonal structure of CeTPO with the space group PJ^/nmm is shown in Fig. [H it 
consists of alternating layers of face sharing TP 4 and OCe4 tetrahedra. The two dimensional 
nature of this lattice in which the CeO layers are well separated by TP layers, together with 
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the presence of conduction electrons, opens up a lot of possibilities for interesting electron- 
mediated magnetic interactions. 

B. CeRuPO, a FM Kondo lattice system 

Magnetic measurements were performed in the temperature range from 2—400 K in a com- 
mercial Quantum Design (QD) magnetic property measurement system (MPMS) equipped 
with an RSO option. The polycrystalline CeRuPO material was powdered, mixed with 
paraffin, and aligned in a magnetic field of ^oH = 5 T at 80°C. NMR measurements dis- 
cussed below suggest that this orientation corresponds to the direction H || ab. In Fig. [2l 
we show the dc-susceptibility x as a function of temperature T at three different magnetic 
fields. The transition into a FM ordered state is clearly visible by a sharp increase of the 
magnetization at Tq = 15 K. The inverse susceptibility shown in the inset of Fig. [2] supports 
that this transition comes from local Ce 3+ moments because the data for T > 100 K can be 
fitted with a Curie Weiss law with an effective moment of fi e g = 2.3 ns/Ce close to that of 
the free ion moment (2.54 fi^/Ce). The Weiss temperature 0\y obtained from this linear fit 
is 0w = +8K. At lower temperatures (T < 50 K), we observe a change of the slope, result- 
ing again in a linear behavior with the fit parameters fJ 1 ^ 7, = 2.2 n^/Ce and 0^r T = +16 K, 
indicating FM correlations. The change of slope and of 0\y can be attributed to the effect 
of the crystal electric field (CEF). A FM state is also evidenced by the magnetization mea- 
surements at 2 K on the same sample, shown in Fig. [31 A small but well defined hysteresis 
can be discerned (see inset of Fig. [3]). However, for the oriented powder the saturation 
magnetization /i sat = 1.2 /is per Ce atom is only reached at 1.2 T. Also the peak observed 
at low field in the susceptibility measurements does not correspond to the expectation for 
a simple ferromagnet. The height of this peak varies from sample to sample but is not 
solely due to domain reorientation because it was also observed in a field cooling experiment 
(open symbols in Fig. [2]). Presently, it is not clear whether these observations are the result 
of measuring powder samples of a system with a complex anisotropic behavior, like e.g., 
a competition between CEF and exchange interaction as in YbNiSn; 4 or if they point to 
a more complex magnetic structure, e.g., a spiral with a very small propagation vector Q. 
Up to now, we could not succeed in growing single crystals large enough to measure the 
orientation dependence of physical properties, to determine the magnetic anisotropy of this 
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system. However, similar behavior of the susceptibility was observed in the FM compound 
TbZrSb, crystallizing in a related structure.— 'i 4 - 

AC-transport resistivity measurements were performed in a standard four-probe geometry 
using a commercial QD physical property measurement system (PPMS). The powdered 
material was pelletized and subsequently annealed at 1200°C to improve crystallinity. The 
temperature dependence of the resistivity is plotted in Fig. HI Three different temperature 
regions can be identified: (I) Above 50 K, the sample shows a linear behavior typical of a 
conventional metal. (II) Below 50 K, there is a pronounced decrease, deviating from the 
linear behavior which is a distinct feature of a Kondo lattice system. This decrease can 
be explained by the onset of coherent Kondo scattering due to the hybridization between 
the localized 4/ and the conduction electrons. (Ill) The FM order is clearly visible in a 
distinct anomaly at Tq = 15 K. In the ordered state, the resistivity follows a power law 
dependence with p oc T 4 . All these features are reproducible, and the very good sample 
quality (residual resistivity po = 1-5 pficm) validates the intrinsic nature of these anomalies. 
In the lower inset of Fig. HI the resistivity of the nonmagnetic reference sample LaRuPO 
is shown which is typical for an usual metal without any anomalies down to 0.5 K. Our 
aim to visualize the high temperature Kondo maximum in the resistivity of CeRuPO by 
subtracting the nonmagnetic data from the resistivity of the CeRuPO sample failed because 
the room temperature resistivity value of LaRuPO is with 700 fiQcm one order of magnitude 
larger than the value of CeRuPO and makes a subtraction very arbitrary. This higher 
resistivity value of LaRuPO is not intrinsic, but can be ascribed to the granularity of the 
sample. However, the presence of Kondo interaction in CeRuPO is further supported by 
thermopower measurements which will be discussed below. 

The specific heat was determined with the PPMS using a standard heat-pulse relaxation 
technique. The 4/ contribution to the specific heat of CeRuPO was obtained by sub- 
tracting the measured specific heat of the nonmagnetic reference sample LaRuPO (C La ) 
from the specific heat of the CeRuPO polycrystals (C Ce ) as C 4 ? = C Ce — C La . The result is 
shown in the upper panel of Fig. [5j plotted as C 4 ^/T vs T for different magnetic fields. A 
large A-type anomaly is evident at Tq = 15 K, indicating a second order transition into the 
magnetic state. The FM nature of the transition is supported by the field dependence of the 
A-type anomaly. The transition temperature shifts to higher temperatures with increasing 
magnetic field; the maximum gets broader and reaches 20 K at HqH = 7T. The extrapolated 
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value of the linear part of the total specific heat C Ce /T at T = K gives 7^ e = 77mJ/molK 2 , 
which is enhanced by a factor of 20 compared to the value 7g a = 3.9 mJ/molK 2 that we 
observe in LaRuPO. This enhancement can be attributed to the 4/-correlation effects. The 
entropy gain at 20 K is i?ln2 and is slightly shifted to higher temperatures if magnetic 
field is applied. At higher temperatures (T > 30 K), does not go to zero, indicating 
further excitations from the 4/ electrons. This contribution to the specific heat is shown 
in the inset and can be explained by a broad Schottky anomaly, resulting from excitations 
to CEF levels of two excited Kramer's doublets at 6 and 30 meV above the doublet ground 
state. For Ce 3+ with J = 5/2 in a tetragonal environment, one expects for the CEF states a 
ground state doublet and two excited doublets. The entropy confirms low lying CEF levels 
because at 150 K already 90% of i?ln6 is reached, which is the total entropy of the whole 
multiplet (see lower inset of Fig. [5]). The absence of a plateau in the entropy around i?ln2 
indicates that the first excited CEF level is not far above Tq. From the reduced jump height 
of the specific heat at the ordering temperature C 4 ^ — C Schottky ~ 9 J/molK, it is possible to 
estimate the ratio of Tk/Tc ~ O.5.— Therefore, the order of magnitude of the Kondo energy 
scale is Tk ~ 10 K, which is in accordance with the onset of coherent Kondo scattering at 
50 K in the resistivity. In the lower part of Fig. [51 we show the temperature dependence of 
the derivative of the resistivity, dp/dT vs T. The anomaly in the resistivity shown in Fig. H] 
is clearly visible in this plot, showing a sharp peak at Tq. This peak behaves very similar 
to that in the specific heat, both at H = and H > 0. Such a pronounced similarity has 
been reported and analyzed by Campoy et a/.— 

To further study the influence of the Kondo interaction on the transport properties of 
CeRuPO, we performed thermal transport (TT) measurements on bar shaped polycrystalline 
samples in a four-probe geometry, using the TT option of the PPMS. Among other transport 
measurements, the thermopower is of particular interest because it is extremely sensitive 
to any variation of the density of states at the Fermi level and exhibits large values (10- 
100 /iV/K) for Kondo lattice systems.— The temperature dependence of the thermopower 
S(T) of CeRuPO is shown in Fig. [6J plotted on a logarithmic temperature scale. S(T) has a 
broad minimum at 150 K (which is negative), followed by a strong increase below 100 K to a 
maximum at 35 K with absolute values as high as 18 /zV/K. A pronounced minimum at Tq 
separate this first maximum from a second one at 10 K. At lower temperatures (T < 10 K), 
the thermopower goes within the error bars linearly to zero. This behavior of S(T) was 
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reproduced in a different sample. Large positive maxima are common in many Ce-based 
Kondo lattice systems and can be ascribed to the Kondo scattering of the different CEF 
levels.— 1 ^ Since the CEF maximum in the thermopower is usually observed around half of 
the splitting energy, 2 - the peak at 35 K in CeRuPO can therefore be nicely attributed to the 
first excited CEF level at 6meV, visible also in the magnetic part of the specific heat. The 
peak position is unaffected by a magnetic field, which confirms the CEF origin. At lower 
temperatures (T < 30 K), the T-dependence of S(T) results from an interplay between 
Kondo effect and FM ordering. Because of the lack of FM Kondo systems, the effect of 
this interplay on S(T) has not yet been studied. The only FM Kondo lattice system where 
thermopower measurements were performed is the related compound CeAgSb2- 21 However, 
CeAgSb 2 presents a strong decrease of S(T) below T c , in contrast to the increase observed 
in CeRuPO. The shift of the minimum in S(T) of CeRuPO to higher temperatures with 
increasing magnetic field indicate that it arises from FM correlations. A likely explanation 
of the temperature dependence of S(T) is that in the absence of FM correlations, S(T) 
would present a broad maximum due to both Kondo effect and CEF splitting,—" 2 ^ as the 
involved energy scales (Tk ~ Tcef) are close to each other. FM inter-site fluctuations lead 
to a reduction of S(T) at Tc, this reduction becoming weaker and shifting to higher T with 
increasing magnetic field. Then, the field dependence of the low temperature maximum 
would just be a consequence of the field dependence of the minimum. To understand the 
observed behavior in more detail, theoretical predictions for S(T, H) in FM Kondo systems 
are necessary. 

The dc-susceptibility x{T) under pressure up to 2GPa was measured in a QD-MPMS 
using a diamond (0.25 carat) anvil cell with a culet size of 1.5 mm. As gasket we used a 
CuBe disk with a bore of 800 //m in diameter. Daphne oil served as pressure transmitting 
medium. The pressure inside the cell was determined by the inductively measured shift 
of the superconducting transition temperature of lead (Pb) placed next to the CeRuPO 
sample, relative to a Pb sample outside of the cell. In Fig. x{T) is shown as function of T 
for different pressures. For 1.4 GPa and 2.0 GPa no significant changes in the temperature 
dependence of x(T) are observed; and therefore, T c stays constant up to 2.0 GPa. In Ce 
based Kondo lattice systems, pressure is expected to suppress magnetic order because the 
volume of the nonmagnetic Ce 4+ ion is smaller. However, it was observed that in a FM 
Kondo lattice system there can be a broad pressure range where Tc only slightly changes.— 
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Further experiments with an extended pressure range are needed to study the influence of 
pressure in CeRuPO in more detail. 

C. CeOsPO, a AFM Kondo lattice system 

Another way to suppress magnetic order is to dope the system chemically with an iso- 
electronic element which leads to stronger hybridization and favors the Kondo interaction. 
Since 5d metals usually cause stronger hybridization than Ad metals, we have synthesized 
CeOsPO. As discussed above, the volume of the unit cell changes only slightly and is 0.5% 
larger than the unit cell of CeRuPO. The inverse susceptibility of CeOsPO is shown in the 
upper inset of Fig. [HI The Curie Weiss fit with /i e fj- = 2.45 /ie/Ce close to that of the free ion 
moment for T > 100 K again shows that Ce is in the trivalent state. The Weiss temperature 
0\v = —9 K obtained from this linear fit indicates AFM correlations. The temperature 
dependence of x (Fig- E]) is typical for an AFM ordered system because the magnitude of 
the susceptibility is nearly field independent. Further on, a peak marks the onset of AFM 
order at Tn = 4.4 K for fi H = 0.2 T. For higher magnetic fields, this ordering tempera- 
ture decreases; for clarity, we only show the data points for three selected magnetic fields. 
At T = 2 K the magnetization increases linearly with the magnetic field, no metamagnetic 
transition is visible up to 5 T (see lower inset of Fig. [S]). 

In Fig. El we present the specific heat C of CeOsPO, plotted as function of temperature. 
We did not subtract any phonon contribution which at low temperatures is negligible com- 
pared to the large mean field type anomaly at = 4.4 K, indicating that the magnetic 
order is a bulk property of CeOsPO. The entropy gain involved in the magnetic transition 
is close to Rln2 at 10 K. The transition becomes broader if magnetic field is applied and 
Tn shifts to lower temperatures, in accordance with the AFM nature of the transition. No 
other transition has been found up to 150 K as can be seen in the inset of Fig. [9] which shows 
C in the whole measured temperature range. Up to now, we were not able to perform any 
reproducible transport measurements of CeOsPO due to difficulties in preparing a compact 
sample of sufficient size. 
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D. 31 P NMR on CeTPO 

We confirmed the conclusions drawn from thermodynamic and transport measurements, 
using results obtained with a microscopic technique, namely 31 P nuclear magnetic resonance 
(NMR). These measurements were performed on polycrystalline samples with a standard 
pulsed NMR spectrometer at a fixed frequency of about 70 MHz and in the temperature 
range between 2 and 300 K. The NMR field sweep spectra were obtained by tracing the 
intensity of the echo as a function of field, using a common spin-echo sequence. NMR 
shift values were determined by using H3PO4 as a reference (K = 0, line in Fig. [TUl) . The 
measurement of the spin-lattice relaxation time was carried out with a saturation recovery 
sequence. The 31 P spectrum (I = |, 7 = 17.1027 MHz/T) in Fig. HDj shows a typical NMR 
powder pattern expected for axial symmetry (tetragonal in this case). The simulation of 
the spectra at a fixed temperature allows the determination of the anisotropic Knight shift 
components 31 K ab (T) and 3l K c (T), corresponding to the H _L c and H || c directions in the 
crystal. The assignment of these components is clearly given by the fact that the highest 
intensity in the powder pattern is attributed to the H _L c direction. A plot K ab vs \ab 
(measured on oriented powder in a magnetic field of 4T), see inset of Fig. [101 evidences a 
linear dependence, as expected from the relation K ab = (A h{ / ^b^ v a)Xo6-~ From the slope 
we determined the hyperfine coupling constant 31 A hi for CeRuPO, 31 A M = 6kOe//iB; a 
value comparable to 31 A hi = 8.73kOe//iB reported for CeP another P containing Kondo 
system.— 31 K a b(T) is positive and reaches its maximum value of about +17% at 2K. When 
the FM order sets in (T C 4T ~ 18 K), the line broadens, but the mean value of the shift stays 
constant with decreasing temperature which is expected for a ferromagnet. In contrast to 
CeRuPO, 31 K ab for the Os compound is negative, and a hyperfine coupling constant of 
31 A^ = — 2.4kOe//iB could be determined— ^ 

Spin-lattice relaxation measurements were carried out at the 31 K ab (T) position in the 
spectra (see Fig. [TUl) . The spin relaxation rate T 1 _1 (T) shown in Fig. (TT] is quite similar 
in the two compounds at high temperatures, while below 40 K, it differs significantly in 
accordance to the different behavior observed in macroscopic measurements. At 300 K, 
Ti has almost the same value in both systems. With decreasing temperature, Tf first 
increases slightly with a smaller slope in CeRuPO than in CeOsPO where it approximately 
follows a power law T^ 1 oc T _1 . Such a behavior is typical for well localized 4/ systems for 
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T » T K && In CeRuPO, Tf l (T) presents a well defined maximum around 35 K before 
decreasing significantly at lower temperatures and merging in a linear in T behavior below 
10 K. Since the temperature of the maximum corresponds to the onset of coherent Kondo 
scattering in p(T), this maximum could be attributed to the onset of the Kondo effect. 
However, the observation that this maximum shifts to lower temperature and becomes more 
pronounced with decreasing magnetic field indicates that it is also related to the onset of FM 
correlations.— A separation of both effects is presently not possible. In contrast, in CeOsPO 
Tf 1 becomes T-independent between 60 and 6 K before it starts to decrease towards lower 
temperatures. Such a plateau like behavior is expected in Kondo lattice systems in the 
temperature range just above Tk-— & This would suggest a Kondo scale of the order of 5 K 
in CeOsPO. The decrease of Tf 1 in CeOsPO below 6K ends in a pronounced drop at the 
temperature corresponding to Tn at the field strength of the NMR experiment. The absence 
of a critical fluctuation induced maximum in Tf 1 at T c in CeRuPO and at T N in CeOsPO 
might be due to the broadening effect of the applied field in the former one and/or to the 
cancellation of these fluctuations at the P-site in the latter one. We included in Fig. [11] 
the temperature dependence of Tf 1 of the nonmagnetic reference compound LaRuPO. A 
common approach in the analysis of T\ in rare earth transition metal compounds is to 
separate the contributions to the spin relaxation, in a contribution of the local moment 
fluctuations which is related to the bulk susceptibility, and a contribution of the itinerant 
conduction electrons which follows a Korringa behavior Tf oc T. This approach seems 
to be quite appropriate for the present compounds, since the difference between Tf *(T) of 
CeRuPO and LaRuPO first increases strongly with decreasing temperature passes through 
a maximum around 30 K before vanishing at low T. A more precise analysis of the NMR 
results shall be performed in a paper devoted to NMR experiments.^ 

III. BAND STRUCTURE CALCULATIONS 

To gain deeper insight into the electronic structure of CeTPO, band structure calculations 
were performed using the full-potential local-orbital minimum basis code FPLO (version 
5.00-19)^ within the local (spin) density approximation (L(S)DA). In the scalar relativistic 
calculations the exchange and correlation potential of Perdew and Wang^ was employed. 
As basis set, Ce (4p4d4f5s5p / 6s6p5d), Ru (4s4p / 5s5p4d), Os (4p4d4f5s5p / Qs6p5d), P 
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(2s2p / 3s3p3d), and O (2s2p3d) where chosen as semi-core/valence states. All lower lying 
states were treated fully relativistic as core states. The inclusion of semi-core states was 
forced by their non-negligible overlap due to the large extension of their wave functions. 
The formally unoccupied P and O 3d states were included to improve the completeness of 
the basis set. The extension of the valence orbitals is controlled by an additional confining 
potential (r/r ) and gets optimized with respect to the total energy— To treat the strong 
Coulomb correlation of Ce 4/ electrons in mean field approximation, LSDA+[7 (in the 
around mean field double counting scheme) was applied with U = 6.1 eV and J = 0.7 eV 
(F 2 = 8.34eV, F 4 = 5.57eV and F 6 = 4.12eV).— In the ferromagnetic cell 546 fc-points, 
and in the antiferromagnetic a/2 ■ a/2 super cells 640 /c-points in the irreducible part of 
the Brillouin zone (BZ), respectively, were used as converged fc-mesh. These /c-mesh sizes 
ensures a numerical accuracy of the total energy at least one order of magnitude better than 
change in energy for the investigated physical quantities. Because of the very small energy 
scales between different magnetic states, CeRuPO was also calculated with the LAPW 
code WIEN2k— for comparison with an alternative full potential scheme.— The resulting 
electronic densities of states, band structures, and total energy differences were basically 
identical for these two band structure codes. 

For both the Ru and the Os compound, we calculated the total energies and the elec- 
tronic structures for three different magnetic states: (I) fully ferromagnetic arranged Ce 
atoms (FM), (II) ferromagnetic arranged Ce atoms in the plane but antiferromagnetic be- 
tween the planes (A-AFM), and (III) antiferromagnetic arrangement of the Ce atoms in 
the plane and between the planes (G-AFM). For CeRuPO, we found the FM state to be 
lowest in energy, followed by the G-AFM (+94.6 me V) and the A-AFM state (+94.9 meV). 
In contrast, for CeOsPO the G-AFM state was slightly favored in energy, followed by the 
A-AFM (+0.08 me V) and the FM state (+1.2 meV). This result is in nice agreement with 
the experimental observation although the rather small energy differences for the Os com- 
pound should be taken with care, especially for the difference between the two AFM order 
types. The corresponding densities of states (DOS) and band structures (FM for CeRuPO 
and G-AFM for CeOsPO) are shown in Fig.[U 

At first glance both band structures and DOS's look very similar. Note that for the Os 
compound (see Fig. [121 upper panel) the number of bands is doubled due to the doubled 
number of atoms in the magnetic super cell, resulting in a folding of the band structure (r 
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and M are identical in the chosen notation of symmetry points referring to the original non- 
magnetic unit cell). The filled valence states originate essentially from P and O p and from 
the transition metal states. The states between about -2eV and 0.5 eV are dominated by 
the transition metal d electrons. In the LDA calculations (not shown), the strong Coulomb 
repulsion of the Ce 4/ electrons is largely underestimated, resulting in an unrealistic peak 
of the partial Ce 4/ DOS at the Fermi level. The application of LSDA+[7 splits the Ce 
4/ states and yields the (qualitatively) correct physical picture of a Ce 3+ state in both 
compounds where the occupied and the unoccupied Ce 4/ states show a split of the order of 
U. Thus, the magnetic Ce spin moment is about 1 hb for both compounds. All other atoms 
show negligible spin polarization, including Ru and Os. This can be seen from the very small 
split of the Ru related bands between -2eV and 0.5 eV (see Fig.[T21 lower panel), especially 
close to the Fermi level. This picture is basically independent from the actual choice of the 
parameter U (in a reasonable range). Fig.[T3l shows the partial 4d DOS for Ru for three 
representative U values (6.1eV±leV).— The dependence on the chosen U is negligible 
in the region around the Fermi level relevant for the low lying excitations, especially the 
metallic properties. This justifies the application of the LSDA+[7 approximation although 
the actual value of U may vary slightly. Therefore, in a first approximation the magnetic 
interaction in the two compounds can be understood from an RKKY type mechanism with 
rather localized Ce 4/ states coupled by the itinerant transition metal d electrons. On the 
other hand one should keep in mind that all dynamic correlation effects are neglected in 
the LSDA+[7 approximation and would require a more sophisticated theoretical treatment. 
Fig. HH shows the detailed band structure close to the Fermi level for CeRuPO (lower panel) 
and CeOsPO (upper panel) together with the transition metal d band characters. Note 
again the doubled number of bands for the Os compound due to the magnetic super cell, 
otherwise the resulting bands are very similar. Neglecting the small spin split for the Ru 4d 
related bands in CeRuPO, we observe a very close similarity with the related superconductor 
LaFePO. 11 

The characteristic feature in the CeRuPO and CeOsPO band structures is the strongly 
pronounced two dimensionality. The in-plane dispersion (along T-X-M-r) exceeds the out 
of plane dispersion (along T-Z) by far more than an order of magnitude. The underlying 
reason is the layer type crystal structure, resulting in a rather weak hybridization between 
adjacent metallic T-P layers. The interaction along the z direction is mainly mediated by 
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the O 2p orbitals of the otherwise non-metallic Ce-0 layers. On the other hand, the slightly 
different oxygen hybridization in both compounds is most likely the reason for the different 
magnetic ground states. Whereas in CeOsPO the bands crossing the Fermi level are mostly 
dispersionless along T-Z, one of the bands in CeRuPO shows sizable dispersion along the 
z direction (see Fig. [H|) . The corresponding band in the Os compound is slightly lower in 
energy and therefore less relevant for the magnetic interaction. Compared to the other T- 
derived bands, the T-d character of this band is somewhat reduced due to the hybridization 
with P and O p electrons. Due to the two dimensional electronic structure, the resulting 
Fermi surfaces (not shown) are almost cylindrical tubes in the case of CeOsPO. In CeRuPO 
one of the cylinders closes due to the dispersion along the z-direction and forms an ellipsoid 
around the Z point like in LaFePO.— To understand the magnetic properties of the two 
compounds in more detail, further more sophisticated studies, including the investigation of 
the role of nesting properties, are necessary. 



IV. DISCUSSION 



Our results on CeRuPO and CeOsPO evidence three interesting features which shall 
now be discussed: The FM ordering and the Kondo behavior in CeRuPO, and the quasi two 
dimensionality. In the absence of more detailed information because of lack of single crystals 
or larger samples for neutron scattering experiments, a deeper insight into the magnetic 
exchange interaction in these systems can be gained by a comparison with RTX compounds 
(R = rare earth, X = Si, Ge) crystallizing in the CeFeSi structure. The arrangement of the R, 
T, and X atoms in this structure is identical to the presented compounds, the only difference 
is the insertion of O atoms between adjacent Ce-planes. B. Chevalier and B. Malaman 36 
discussed the magnetic order in these RTX compounds using a simple model with three 
exchange interactions: Jo within an R-plane (this corresponds to the next nearest R-R- 
neighbor), Ji between R atoms separated by a TX-layer, and J2 between adjacent R-planes 
sandwiching an O-layer (this corresponds to the nearest R-R- neighbor). Systematic studies 
of these RTX compounds showed that Jo is always positive (FM), while Ji and J2 change sign 
depending on the T element. J\ and J2 are positive in the RFeSi light rare earth compounds, 
while Ji < 0, Ji > for T = Co, Ti, and Ji > 0, J2 < for T = Ru. The only exception 
with both Ji and J2 negative is CeCoGe. Thus, for this type of arrangement of R, T and 
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X atoms there is a clear predominance of FM exchange. Notabene strong ferromagnetism, 
with Tc = 374 K in GdTiGe, has also been observed in the related CeScSi structure which 
differs only in the long range stacking of the layers.— Therefore, the observation of a FM 
state in CeRuPO is not surprising and can be related to an atomic arrangement promoting 
FM interactions. Comparing the FM state in CeRuPO with the AFM structure found in 
RRuSi(Ge), the inclusion of an O atom between adjacent Ce-layers switches J2 from AFM to 
FM. In absence of more detailed information, it is yet not clear whether the AFM structure 
in CeOsPO is due to a negative J\ or a negative J2. A preliminary analysis of the change in 
the anisotropy of the hyperfme coupling constant between CeRuPO and CeOsPO suggest 
that at least J\ is becoming AFM in the latter one since P is along the corresponding 
exchange path.— 

As already stated, only very few FM Kondo lattices have been reported previously. None 
of them show such a strong decrease in the resistivity starting far above Tc as we observed in 
CeRuPO. As an example YbNiSn, CeAgSb 2 ,— and CeSii.71— present a weak increase in p(T) 
below 50 — 100 K ending in a maximum at a temperature T max which is at or only slightly 
above Tc. Thus p{Tc) is om y slightly smaller than p(T max ). In contrast, p{T) in CeRuPO has 
at Tc decreased to less than 50% of the value at 50 K which would roughly correspond to the 
position of the maximum after subtraction of the phonon scattering. Thus spin flip scattering 
is already significantly suppressed at Tc, indicating that the Kondo fluctuations at Tc are 
likely much stronger than in the previously known FM Kondo lattices. The most appropriate 
compound for a detailed comparison is the closely related CeAgSb2 38 which crystallizes in 
the same structure type as CeTPO. However, the much larger size of Sb compared to O 
and P leads to quite different Ce-Ce, Ce-T, and Ce-P distances. Nevertheless, CeAgSb 2 
also orders ferromagnetically at Tc = 9.6 K. While at ambient pressure p(T) in CeAgSb2 
decreases by only a few percent between T max and Tc under an isostatic pressure larger than 
2.4 GPa, its temperature dependence becomes comparable to that in CeRuPO at p = 0; 
although, the onset of coherence stays at a much smaller temperature, indicating a smaller 
Tk. The ratio RRt c = p(T max ) / p(Tc) = 2.2 in CeAgSb 2 at p = 3.26 GPa which is the highest 
pressure where the FM transition is still clearly visible is close to RRt c = 2.3 in CeRuPO at 
ambient pressure, suggesting a comparable T K /T C ratio under these conditions. However, 
because of the much lower Tc = 2.4 K at p — 3.26 GPa in CeAgSb2 this also corresponds 
to a much lower absolute value of Tk- A smaller Tk in CeAgSb2 is further evidenced by 
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the weaker thermopower.— The larger Tk in CeRuPO is in accordance with general rules 
for the strength of hybridization between / and conduction electrons first proposed by D.D. 
Koelling et al.— Both the smaller size of P compared to that of isoelectronic Sb, and the 
shift of the ci-states closer to the Fermi level upon exchanging Ag by Ru are expected to 
lead to a stronger hybridization. One can also suspect that the replacement of one Sb by 
O has the tendency to push Ce towards a higher oxidation state. The effect of pressure on 
CeRuPO and the way in which a further increase of Tk suppress the FM order shall be an 
important issue for future experiments. 

Quasi two dimensional (2D) intermetallic Kondo lattices have attracted considerable at- 
tention in the past years. The main reason was the discovery of a rather high superconducting 
transition temperature T c = 2.3 K in CeCoIn 5 ; the layered structure of this compound is 
suspected to be at the origin of the enhanced T c .— However, the 2D character of these 115 
compounds is not very pronounced; the slope of the upper critical field at T c e.g. shows only 
a rather small anisotropy of a factor of two. Only one of the three sheets of the Fermi surface 
of CeCoIns has a clear 2D character while the two others are clearly 3D. 42 Thus, intermetal- 
lic heavy Fermion compounds with a better defined 2D character are of strong interest. The 
results of our LDA+[7 calculations as well as those on LaFePO^ suggest that the electronic 
states at the Fermi level of the presented RTPO have a stronger 2D character than in the 
CeTIn 5 compounds, since at least four of the five sheets of the Fermi surface are 2D. Thus 
the growth of single crystals large enough for allowing the study of the anisotropy of the 
transport and electronic properties shall be another important issue for future research. 



V. CONCLUSIONS 



In summary, we have shown that the CeTPO (T = Ru, Os) compound series presents an 
interesting class of Ce-based Kondo lattice systems at the border between intermetallic and 
oxide materials with a layered structure and strong tendency to ferromagnetism. CeRuPO 
seems to be a nice example of a FM Kondo lattice (Tc = 15 K, Tk ~ 10 K). Together with 
the antiferromagnetic CeOsPO (T$ = 4.4 K), these two compounds are ideally suited to 
study the difference between FM and AFM Kondo lattices. This change from the FM state 
in CeRuPO to the AFM state in CeOsPO is induced by only minor changes in the crys- 
tallographic parameters. However, further pressure experiments on CeRuPO are necessary 
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to investigate the interplay between Kondo interaction and development of long range mag- 
netic order. The yet unique combination of quasi two dimensional electronic properties and 
strong FM correlations could lead to new kind of quantum critical behavior. Furthermore, 
a systematic study of related derivatives by replacing the T atoms with other d elements 
brings up a lot of possibilities for further interesting materials. 
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FIG. 1: (Color online) Tetragonal crystal structure (P4/nmm) of the CeTPO compound series, 
showing the alternating layers of TP4 and OCe4 tetrahedra. 
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FIG. 2: (Color online) x vs T a t various magnetic fields along the easy direction of oriented powder 
of CeRuPO. The magnetic transition into a FM ordered state is clearly visible at Tq = 15 K. 
Open symbols indicate data from a field cooled experiment. The inset shows the inverse magnetic 
susceptibility together with a Curie Weiss fit (/x e ff = 2.3 /xb and ©w = +8K). 
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FIG. 3: Isothermal magnetization as a function of applied magnetic field at T = 2K. The inset 
shows the small hysteresis at smaller magnetic fields. 
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FIG. 4: Temperature dependence of the resistivity of CeRuPO. The pronounced decrease below 
50 K is attributed to coherent Kondo scattering (see upper inset) . The magnetic order is visible as 
an additional decrease at 15 K. In the lower inset, the temperature dependence of the resistivity of 
the nonmagnetic LaRuPO is shown. 
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FIG. 5: (Color online) (a) Temperature dependence of the magnetic part of the specific heat divided 
by temperature C 4 ^/T of CeRuPO. The sharp A-type anomaly at Tq = 15 K corresponds to the 
magnetic order and shifts to higher temperatures in applied magnetic field. The inset shows C 4 ^ 
at higher temperatures, the broad Schottky peak above the transition is due to excitations of the 
crystal electric field levels, (b) Derivative of the resistivity of CeRuPO with respect to temperature 
plotted as dp/dT vs T, being very similar to the specific heat data C^/T. The inset shows the 
temperature dependence of the entropy, calculated from specific heat data. 
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FIG. 6: (Color online) Thermoelectric power vs temperature of CeRuPO for different magnetic 
fields. 
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FIG. 7: (Color online) \ vs T measured in a magnetic field of 1 T at atmospheric pressure and 
under an applied pressure of p = 1.4 GPa and 2 GPa, respectively. The data for the different 
pressures are shifted with respect to each other for clarity. 
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FIG. 8: (Color online) x vs T f° r three different magnetic fields along the easy direction of oriented 
powder of CeOsPO. The magnetic transition into a AFM ordered state is visible at Tn = 4.5 K. The 
upper inset shows the inverse magnetic susceptibility together with a Curie Weiss fit (/x e ff = 2.45 
and 0w = — 9K). In the lower inset, the magnetization is plotted vs magnetic field at T = 2K. 
No metamagnetic transition is visible up to 5 T. 
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FIG. 9: (Color online) Temperature dependence of the specific heat C of CeOsPO. The phonon 
contribution is not subtracted in this representation. The large anomaly at Tn = 4.4 K corresponds 
to the magnetic order and shifts to lower temperatures in applied magnetic field. The inset shows 
the specific heat of CeOsPO in an extended T range. 
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FIG. 10: (Color online) Spectrum of CeRuPO at T = 100 K and v = 70 MHz. The solid line 
describes a simulated spectrum. The inset shows the linear behavior of K a b vs Xab- 
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FIG. 11: (Color online) Spin lattice relaxation rate 31 (1/Ti) of CeRuPO (filled circles) and CeOsPO 
(open circles), together with 31 (1/Ti) of the nonmagnetic LaRuPO (filled squares). 
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FIG. 12: (Color online) Band structure and total as well as partial Ce densities of states for the FM 
state of CeRuPO (lower panel) and the G-AFM state of CeOsPO (upper panel). The symmetry 
points of the BZ correspond to the primitive cell for better comparability. The Ce contributions 
to the densities of states are indicated by different fillings. The spin down bands for CeRuPO are 
indicated by dashed (red) lines with small circles. 
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FIG. 13: (Color online) Partial Ru Ad DOS of CeRuPO for three different choices of the Coulomb 
repulsion parameter U. The inset shows a narrow region around the Fermi level where the curves 
are identical within the line widths. Small differences only occur in energy regions with significant 
Ce 4/ contributions (see Fig. [T2j) . 
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FIG. 14: (Color online) Band structure close to the Fermi level for CeRuPO (lower panel) and 
CeOsPO (upper panel). The symmetry points of the BZ correspond to the original non- magnetic 
cell. The transition metal d band characters are indicated by the size of the filled circles. The two 
spin directions for Ru are shown by black and grey (red) symbols. 
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